Probe spectroscopy in an operating magneto-optical trap: 
the role of Raman transitions between discrete and continuum atomic states 
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We report on cw measurements of probe beam absorption and four-wave-mixing spectra in a *^Rb 
magneto-optical trap taken while the trap is in operation. The trapping beams are used as pump 
light. We concentrate on the central feature of the spectra at small pump-probe detuning and 
attribute its narrow resonant structures to the superposition of Raman transitions between light- 
shifted sublevels of the ground atomic state and to atomic recoil processes. These two contributions 
have different dependencies on trap parameters and we show that the former is inhomogeneously 
broadened. The strong dependence of the spectra on the probe-beam polarization indicates the 
existence of large optical anisotropy of the cold-atom sample, which is attributed to the recoil 
effects. We point out that the recoil-induced resonances can be isolated from other contributions, 
making pump-probe spectroscopy a highly sensitive diagnostic tool for atoms in a working MOT. 

PACS numbers: 32.80.Pj, 42.50.Vk, 42.65.-k 



I. INTRODUCTION 

The magneto-optical trap (MOT) is now a standard 
tool for obtaining cold atomic samples. Due to signif- 
icantly reduced Doppler broadening, low collision rate 
and long interaction time, such samples offer a unique 
possibility of ultra-precise spectroscopic measurements. 
They are especially useful for experiments which inves- 
tigate Raman transitions between nearly degenerate en- 
ergy levels. The resonant structures in probe spectra 
observed in cold media for pump-probe detuning much 
smaller than the natural linewidth are due to numer- 
ous effects: stimulated Rayleigh scattering Q and Ra- 
man transitions between either Hght shifted sublevels of 
a ground atomic state |2], vibrational energy levels of 
atoms localized in an optical lattice 0|, or kinetic mo- 
mentum states of unbound atoms 0|. However, the MOT 
is not an ideal tool for systematic investigation of these 
effects. Because of the fixed three-dimensional geometry 
of the trapping beams, the presence of inhomogeneous 
trapping magnetic field, and limited fiexibility in varying 
the trap parameters, it is difficult to selectively address 
the phenomena mentioned above. This is why the MOT 
usually serves only as an initial stage of preparation of a 
cold sample which, after switching off the MOT's opti- 
cal and magnetic fields, is subject to more precisely con- 
trolled experimental conditions, such as dedicated pump 
beam geometry and polarization phase- 
stabilized pump beams [Sj, etc. 

Despite these difficulties, cw pump-probe spectroscopy 
of trapped atoms performed while the MOT is working, 
using the trapping beams as pump Hght, has been studied 
by several authors. Such experiments were first carried 
out by Tabosa et al. P| and Grison et al. 0. Also, there 
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exists spectroscopic evidence of optical lattices in a MOT 
with a special, phase-shift-insensitive, trapping beam ge- 
ometry [ly|. These experiments indicated the potential of 
probe spectroscopy for cold-atom diagnostics. However, 
the complexity of such spectra from a working MOT re- 
sulted in their possible appHcation for MOT diagnostics 
being neglected. In this paper we revisit the problem and 
describe our systematic studies of such spectra. 

The aim of this paper is to examine the Raman pro- 
cesses which occur in a MOT. Their interplay makes 
the MOT spectra far more complex than those obtained 
under simpler experimental conditions. The systematic 
study of the MOT spectra is necessary for their possible 
application for MOT diagnostics which, as we show be- 
low, can be successfully performed by pump-probe spec- 
troscopy. Working MOT spectroscopy can constitute a 
powerful, non-destructive diagnostic tool. It can pro- 
vide information on the atomic cloud density and 
average Rabi frequency Below, we concentrate on 

the additional possibility of spectroscopic velocimetry of 
atoms in a working MOT, based on the recoil-induced 
resonances (RIR). This method has been already success- 
fully used to determine the temperature of the cold sam- 
ple but under conditions in which RIR 
constituted the only relevant contribution to the spectra. 
In a working MOT, however, recoil-induced resonances 
appear accompanied by the other Raman processes which 
necessitates clear identification and more thorough anal- 
ysis of their individual contributions. 

Below, in Sec.|n]we briefiy characterize our setup and 
in Sec, mil present in more detail our experiment and re- 
sults. We perform spectroscopic measurements of ^^Rb 
atoms in a standard MOT 01 • With the MOT's opti- 
cal and magnetic fields still turned on, we simultaneously 
record two signals: a probe absorption and a four-wave 
mixing signal. The latter is generated by atoms as a 
result of their nonlinear interaction with the probe and 
trapping beams and in such a geometry propagates oppo- 
sitely to the probe beam direction The spectroscopic 
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signals were acquired for various trap parameters, e.g., 
the trapping beams' intensity and detuning from atomic 
resonance and the gradient of the MOT's magnetic field. 
This allowed us to carry out a systematic investigation of 
the phenomena that affect the shapes of the absorption 
and four-wave mixing spectra. 

In Sec. IIVL we show that the theoretical explanation 
of the observed spectra has to include two contributions: 
the Raman resonances between light-shifted Zeeman sub- 
levels of the ®^Rb {F — 3) ground state, or for 
short, the Raman-Zeeman resonances (RZR) B. flSll and 
the recoil-induced resonances (RIR) U, [lai, |2d|. Ad- 
ditionally, in order to reach good agreement between 
the theory and the experimental data, one has to con- 
sider realistic physical conditions in the MOT, namely, 
the light intensity and polarization gradients, the three- 
dimensional geometry of the trapping beams, and the 
inhomogeneity of the quadrupole MOT's magnetic field. 
In Sec. IIVI we present interpretation of the main indi- 
vidual mechanisms that can be recognized in the MOT. 
We also note that the difference between the cr+ and 
(7^ absorption spectra is the signature of the optical 
anisotropy of the cold atomic sample in a MOT and is 
due to recoil resonances induced by the probe and dif- 
ferent trapping beams. Finally, we suggest a method to 
eliminate the RZR contribution, thus opening the pos- 
sibility of real-time, nondestructive MOT temperature 
measurement based on observation of the width of RIR 
contribution [3| and conclude our work in Sec.lvl 



II. EXPERIMENTAL SETUP 

We use a standard magneto-optical trap in 
a stainless-steel vacuum chamber with anti-reflection 
coated windows. Our laser setup consists of four home- 
built diode lasers. One of them, equipped with an ex- 
ternal cavity, is frequency stabilized using saturated ab- 
sorption in a rubidium vapor cell, or alternatively the 
Doppler-free dichroism method j22|. This laser serves 
as the master for injection seeding into the trapping 
and probe lasers. The 780 nm light from a 70-mW 
trapping diode laser is divided into three beams, which 
have Gaussian radius a « 0.4 cm and peak intensity 
Imax = 12 mW/cm^. They are retro- reflected after pass- 
ing the trap cell. The frequency shift of the trapping and 
probe beams is controlled by several acousto-optic mod- 
ulators (AOMs). The typical detuning, A — uj ~ ujq, 
of the trapping beams of frequency uj from the trap- 
ping transition resonance frequency uq is — SF, where 
F = 27r X 5.98 MHz is the transition's natural linewidth. 
Before injection into the probe laser diode, the master 
beam double passes the AOM, which allows us to tune the 
frequency ujpr of the probe laser in the range of ±40 MHz 
around the master laser frequency. In the present exper- 
iment the probe sweep range was ±3 MHz. Since both 
trapping and probe lasers are seeded with the same beam, 
they are phase locked. Thus, when these two lasers work 




Figure 1: The geometry of the experimental setup. 



in the pump-probe conflguration they provide a spectro- 
scopic resolution which allows one to distinguish resonant 
structures with widths of a few kHz, limited by residual 
phase fluctuations. The repumping beam is derived from 
an independent, free running laser diode, and is over- 
lapped with one of the trapping beams. The quadrupole 
magnetic fleld is generated by a pair of anti-Helmholtz 
coils, which produce an axial gradient of 16 G/cm. Stray 
dc magnetic flelds are zeroed by three orthogonal pairs 
of Helmholtz coils. We trap about 10^ atoms in the 
cloud with a Gaussian radius of ct « 0.9 mm. The tem- 
perature of our sample, measured by the time-of-flight 
method [l^l, is about 100 ^K. The geometry of our ex- 
periment is depicted in Fig.Q] The probe beam enters the 
cloud of cold atoms at a small angle — 3.5° with one of 
the trapping beams. The polarization of the probe beam 
is set by the A/4 waveplate placed after a polarizer (P). 
We deflne the polarization of the probe beam with re- 
spect to the nearly co-propagating trapping beam, which 
is cr~ -polarized. We call the probe polarization cr^ (cr"'") 
when it is the same as that of the nearly co-propagating 
(counter-propagating) trapping beam. After traversing 
the cloud, the probe beam is directed onto a photodiode 
which records the absorption spectrum (ABS phd). The 
four-wave mixing beam, generated in the cloud and prop- 
agating oppositely to the probe beam, is reflected by a 
50/50 beamsplitter (BS) onto another photodiode (FWM 
phd) . Both signals are acquired simultaneously and thus 
can be directly compared. The probe beam is shaped to 
have the diameter smaller than the cloud size in order 
to avoid an undesired background of non-absorbed light. 
Typically, the probe-laser power is about 1 /xW, and its 
frequency is swept at a rate of 5 MHz/s. The spectra 
presented in this paper are recorded as a function of the 
pump-probe detuning, S — (jjpr — w, and are averaged over 
20 probe sweeps. 



III. MEASUREMENTS AND RESULTS 

We have performed systematic studies of absorption 
and four-wave mixing spectra of *^Rb atoms in a MOT. 
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Figure 2: Probe absorption spectra (a) and four-wave mixing 
signals (b) for the cr"*" and cr~-polarized probe beam. The 
spectra are recorded as a function of the pump-probe detuning 
5. /max stands for the peak intensity of a single trapping 
beam. The trapping beam detuning from atomic resonance 
is A = — 3r, the magnetic field gradient along coil axis is 
dB/dz = 13 G/cm. 



We have varied the trapping beams' intensity, their de- 
tuning and the magnetic field gradient. Examples of 
absorption and four-wave-mixing spectra simultaneously 
recorded for three various trapping beam intensities are 
depicted in Fig. |2l The shapes of the spectra are sim- 
ilar to those previouslvobserved and assigned to RZR 
(see, for example, 0, However, in contrast to the 

pure RZR case, they differ significantly for cr'^ (left col- 
umn) and (J~ polarization of the probe beam (right col- 
umn). In particular, in absorption, the a"*" and <j~ spec- 
tra have different amplitudes and positions of individual 
resonances. In four- wave mixing, the spectra consist of 
two very distinct contributions, one broad and second 
ultra-narrow. The broad contribution to the wave-mixing 
spectra recorded with the ct"*" probe appears broader that 
that associated with the a~ probe. 

From Fig. [21 it is clearly seen that with the increase 
of the trapping beam intensity, both the absorption and 
four-wave-mixing spectra become wider and their reso- 
nant structures better resolved. This fact can be quali- 
tatively explained by the increase of the splitting of the 
ground-state Zeeman sublevels due to the ac Stark shift 
(light shift). For > F, this shift is proportional to the 
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Figure 3: Probe absorption spectra (a) and four-wave mix- 
ing signal (b), recorded when almost all trapping beam power 
is sent to the pair of trapping beams nearly collinear with 
the probe, as schematically shown in (c). The intensity of 
these beams is Iz = 30 mW/cm^, intensity of beams transver- 
sal to the probe beam is 3 mW/cm^. The detuning of the 
trapping beams is A = — 3r and the magnetic field gradient 
is dB/dz = 13 G/cm. The dispersive shape of the recoil- 
induced resonance shows up in the center of the cr~-polarized 
probe absorption spectrum (marked by grey arrow). 



pump beam intensity {ft is the Rabi frequency associ- 
ated with the pump field). The quantitative description 
of the shape of the spectra, however, should incorporate 
all possible Raman transitions between the pairs of adja- 
cent Zeeman sublevels in the (^"=3) ground state of 
^^Rb (see Fig.^ below) and also include a possible RIR 
contribution. This will be discussed in the next section. 

In the case of the absorption spectra for the a~- 
polarized probe, a weak, ultra-narrow structure appears 
near the pump-probe detuning (5 « 0. In the four- 
wave mixing signal this ultra-narrow resonance is even 
more pronounced and occurs for both polarizations of 
the probe beam. This feature in a MOT spectra has 
not yet been thoroughly examined. We attribute it to 
recoil- induced resonances To verify this hypothe- 
sis we have performed an additional measurement. The 
trapping beam intensities were set in such a way that 
almost all trapping laser power was sent to the pair of 
trapping beams nearly collinear with the probe. In such 
a situation, trapping beams transverse to the probe are 
much attenuated and serve only to sustain the cloud sta- 
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bility. Thus, we approach the one-dimensional pump- 
probe spectroscopy setup with two counter-propagating 
strong pump beams nearly collinear with the probe @, 
while still having cold atomic cloud in a stable MOT. 
Spectra registered under such conditions, presented in 
Fig- El reveal resonant structures around (5 w even bet- 
ter resolved than those of Fig. [21 In particular, a distinct 
dispersion-like resonance develops in the center of the ab- 
sorption spectrum for the tT~-polarized probe, while for 
the cr+ probe polarization it is barely visible. The shape 
of this resonance agrees well with the theory of recoil- 
induced resonances and its width corresponds to the 
actual temperature of the atomic sample, which has been 
measured independently (2^ . 

In order to understand how the measured spectra de- 
pend on realistic MOT conditions we have measured 
the absorption and four-wave mixing spectra for vari- 
ous magnetic-field gradients. Since the size of the cloud 
is proportional to {dB / dz)~^/'^ the value of the 

magnetic field in the peripheries of the cloud scales as 
{dB/dzYf^. For a well-aHgned MOT with the value 
B — Q m the trap center, we observed that for smaller 
gradients the resonances occur at the same frequencies 
as for larger gradients but are better resolved. Thus, due 
to the finite size of the atomic cloud, when modeUing 
the experimental curves one has to take into account the 
inhomogeneous broadening due to the MOT quadrupole 
magnetic field. 



IV. INTERPRETATION 
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Figure 4: (a) Multilevel structure of the 5 {F = 3) 

*^Rb ground state. Zeeman sublevels are perturbed with a tt- 
polarized pump (thick grey arrows), which results in their 
energy shift and symmetric population distribution (align- 
ment), marked by grey circles. The structure is then probed 
by a (T^-polarized beam, (b) Theoretical absorption signal 
compared with the experimental spectrum. The theoretical 
curve was generated for the Rabi frequency Q, — GF, detuning 
A — — Sr and 7i_fc — 0.13F for all i, k. Experimental condi- 
tions: /max = 10.2 mW/cm^ A = -SF, dB/dz = 13 G/cm. 



The shape of the absorption and four-wave mixing 
spectra that we measure in our experiment can be ex- 
plained on the basis of two processes: Raman transi- 
tions between Hght-shifted Zeeman sublevels (RZR) of 
the ground state 0, 113 and recoil- induced resonances 
(RIR) 0. The theoretical outline of both processes is 
presented below. However, in order to perform the com- 
plete modelling of the spectra, we extended the present 
theories by inclusion of two types of inhomogeneous 
broadening: one due to modulation of the net E-field 
in the trap and the other connected with the presence of 
the quadrupole trapping magnetic field B. Moreover, the 
standard three-dimensional geometry of the MOT beams 
has to be taken into account when considering recoil pro- 
cesses in a working trap. 



A. Raman transitions between light-shifted 
Zeeman sublevels (RZR) 

Let us consider the cycling 5 ^Si/2 {Fg = 3) — 5 
{Fe = 4) transition in *^Rb atom in a MOT. The atom is 
subject to three pairs of counter-propagating laser beams 
of frequency to with pairwise orthogonal circular polar- 
izations. Their interference in the intersection region re- 
sults in a complex spatial modulation of light intensity 



and polarization j23|. For the sake of simpHcity, we con- 
sider here only the case where the net light field is lin- 
early polarized and choose quantization axis parallel to 
the local direction of the net field E. In this reference 
frame, the resulting pump fight is 7r-polarized and shifts 
the Zeeman sublevels of both the ground and the ex- 
cited atomic state. The 5 ^Si/2 {Fg = 3) ground state 
multilevel structure is presented in Fig. The opti- 
cal pumping in such a scheme leads to alignment, i.e. 
a symmetric distribution of populations with respect to 
the TOg = sublevel of the ground state, with this sub- 
level being mostly populated. The atoms interacting with 
the pump li ght are probed by circularly polarized weak 
probe laser [2a|- The Raman processes involving a tt- 
polarized pump and a cr^ -polarized probe photon lead 
to transitions with Am^ = ±1. For a given population 
distribution among the nig sublevels, two directions of 
such processes are possible, depending on the sign of the 
probe-pump detuning 5. In the considered case of dom- 
inant population in the rUg = sublevel, the Raman 
transitions with Am = — 1 take place for 5 <Q and lead 
to gain of the probe; those with Am — +\ take place for 
(5 > and result in its attenuation. These processes are 
resonant whenever \5\ coincides with the energy separa- 
tion of the adjacent sublevels. The amplitude of the cor- 
responding resonance is proportional to the population 
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difference of the sublevels involved in the transition and 
depends on the Clebsch-Gordan coefficients associated 
with the specific transition path. The simplest model of 
the probe absorption spectrum is obtained by summing 
the Lorentzian profiles centered at the appropriate reso- 
nance frequencies and weighted by products of the rele- 
vant population differences and squares of the Clebsch- 
Gordan coefficients 




srzr(^) = X! Wi-i,iAUi^i^iL {S, (5j_i,i, 7i-i,i) 

i=l 

-3 

- ^ w,+i,iAni+i,iL((5, (5i+i,i,7i+i,.,) . (1) 



i=-l 



In the above equation, Wkj is the weight associated with 
the Clebsch-Gordan coefficient along the k — irig 
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Raman transition path, AUkj is the population dif 



ference of the k 



and j 



sublevels, and the 



Lorentzian profile is given by 



Lis, 6k,j,"fk,j) 



(2) 



where Skj and jkj are, respectively, the resonance fre- 



quencies and widths of the k 



j = nig transitions. 



These widths are due to the optical pumping, the finite 
interaction time and atomic collisions in a trap. The 
spectrum generated using formula 10 with assumption 
of full optical pumping in the closed system shown in 
Fig. 3t is presented in Fig. ^ compared to the relevant 
experimental data. The comparison shows that the lat- 
ter constitutes an envelope for the theoretical curve un- 
der which individual resonances are unrealistically well 
resolved. This suggests a broadening mechanism, which 
is indeed provided by the spatial inhomogeneity of the 
light intensity in the trap. To account for this effect, we 
averaged the calculated spectra over Rabi frequency fl. 
This can be done as described in the Ref. |2ZI]; but in 
our case it was sufficient to sum the absorption spectra 
generated according to Eq. for the appropriate range 
of Rabi frequencies. By numerical simulations of var- 
ious field contributions, similar to that of Ref. j2^, we 
have found that the relevant range that produces the best 
agreement between the theory and the experiment is 5T 
to 7.5r. This corresponds to the average Rabi frequency 
of 6.25r determined as in Ref. The inclusion of the 
E-field inhomogeneity yields a better agreement of theo- 
retical results with the experimental data, as presented in 
Fig. El However, there is still a discrepancy between the 
theoretical prediction and the experimental signal near 
i5 ~ 0. Moreover, the model discussed above does not 
explain the profound difference between the spectra reg- 
istered for tJ"'" and a~- polarized probe. The explanation 
of this effect is presented in the next subsection. 



Figure 5: The comparison of experimental absorption spec- 
trum (acquired under the same experimental conditions as in 
Fig.Et) with the theoretical curve generated according to JJl 
and averaged over the Rabi frequency range SI G [SF; 7.5r]. 



B. Recoil-induced resonances (RIR) 

The momentum exchange between non-localized atoms 
and the laser field is associated with the phenomenon of 
the so-called recoil-induced resonances (RIR) . They were 
predicted by Guo et al. ^2Qj and observed by Grynberg 
et al. Q, 0- A simple momentum-space analysis of RIR 
can be found in Refs. Q, Here we briefly recall 

its results. 

Let us consider atoms of mass m interacting with a 
weak probe beam (with wave vector kpr and frequency 
Lo + S), which makes an angle 9 with the direction of prop- 
agation of a strong pump beam (wave vector k, frequency 
uj), as depicted in Fig.^K. The process of pump photon 
absorption followed by the probe photon emission results 
in the momentum change Ap of an atom, where 



Ap 



-2hk sin — . 
2 



(3) 



The preceding equation is derived assuming small probe- 
pump detuning |k| « |kpr| — k. The resonance occurs 
whenever the probe-pump detuning coincides with the 
kinetic energy difference, namely at detuning 



2k 



hk sin V \ sm ■ . 

2^/2 



(4) 



An analogous consideration performed for the case of 
probe-photon absorption followed by emission of a pho- 
ton into the pump beam leads to the nearly identical for- 
mula for (5res as Q, but with the opposite sign. Similarly 
to the RZR case, the resonance amplitude is proportional 
to the population difference between the relevant atomic 
states, but now we consider the continuum, kinetic mo- 
mentum states rather than discrete, magnetic sublevels 
(Fig.Eb), 



AH (pflnal,Pinitial) = H (p + Ap) - H (p) 



(5) 



Integration over all possible momentum values leads to 
the formula for the RIR signal given as 



srir(<5, ^) 



dpAn(pflnal,Pinitial)-^('5)- (6) 
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Figure 6: (a) Basic geometry for the recoil induced resonance, 
associated with atomic interaction with two beams: the pump 
(thick grey arrow) and probe (thin arrow) beams intersecting 
in an atomic sample at a small angle 9. (b) Raman transitions 
between the kinetic momentum states for an atom in a state 
with specific mg. 



Assuming that n(p) is a Gaussian distribution, Ap <C 
PT = \fmkBT and 7 <C kpT/m 0| (all these conditions 
are fulfilled in our case), one arrives at an analytic for- 
mula for the absorption spectrum signal 



srir((5, 0) 



hS 



X exp 



2uy^ksm (61/2) 
52 



2w2,(2fcsin(0/2))2 



(7) 



where ut = pr/m is the most probable atomic speed. 
The signal Q is the derivative of a Gaussian function. 
It has its minimum/maximum for S — ±2ut^ sin(0/2) 
and the width of the spectrum A^/fl, defined as the dis- 
tance between minimum and maximum, is proportional 
to VT. The RIR spectrum can thus serve as a spectro- 
scopic tool for the ternperature measurement of a cold 
atomic sample 



according to the formula 



T - 



le/csF sin(6i/2) 



^RIR- 



(8) 



The results presented above refer to a two-level atom and 
a simple pump-probe beam configuration. Nevertheless, 
this simple approach can be also successfully used in our 
case in which multilevel atoms are subject to a three di- 
mensional set of trapping beams. Since a distinct RIR 
contribution appears only near S ^ (see Figs. El EJ, 
it must result from Raman transitions between the mo- 
mentum states of the same magnetic sublevel, Airig = 0. 



In principle, RIRs can also occur between momentum 
states associated with various Zeeman sublevels, i.e. with 
Anig 7^ 0, in which case some ja^ asymmetry might 
be visible, but this contribution to the overall signal is 
negligible. First, the amplitudes of such contributions 
have been shown to be about order of magnitude weaker 



than of the Am^ = principal ones 



Second, they 



would result only in a small frequency shift of the broader 
RZR contribution. Third, as shown below in IIV (;l the 
Raman transitions with Am^ 7^ in a working MOT 
are subject to inhomogeneous broadening which further 
reduces their importance for RIR. Thus, the recoil con- 
tribution is due to the probe and pump beams of the 
same polarization. This implies that our atoms can be 
effectively treated as two level systems [1^ and that 
the polarization of the probe selects the pump of the ap- 
propriate polarization among all six available beams. 

For example, consider the case in which the probe 
beam is cr~-polarized and the quantization axis is par- 
allel to its direction of propagation. Obviously, one con- 
tribution to the recoil process is due to the probe beam 
combined with the nearly co-propagating trapping beam. 
Analogously, the (T+-polarized probe interacts with the 
nearly counter-propagating trapping beam. The second 
contribution is due to the trapping beams transverse to 
the probe. In the chosen reference frame, these beams ap- 
pear as linearly a and 7r-polarized. Depending on their 
relative phases, they represent cr+ or cr^-polarized pump 
photons with the same probability, so they equally con- 
tribute to the recoil processes both for the and a~ - 
polarized probe. This consideration yields the total RIR 
signal in the probe beam absorption in the MOT with all 
six trap beams as equal to 



srir, mot('5) 



srir(,5,0) + 2srir((5,9O°) 
for cr~ -polarized probe; 

srir(,5,18O°-0) + 2srir((5,9O°) 
for fj^-polarized probe. 



(9) 



The above discussion allows one to associate the observed 
distinct difiPerence between the MOT spectra taken for 
and fj "'"-polarized probe exclusively with the recoil pro- 



cesses, since, as it was pointed out in Sec. IIV Al the RZR 
contributions are insensitive to the probe beam polariza- 
tion. 



C. Influence of the MOT magnetic fleld 

The considerations presented above lead to the conclu- 
sion that both RZR and RIR contribute to the observed 
spectra, so the complete formula for the absorption signal 
becomes 



s((5) = srzr(^) H- srir,mot(^), 



(10) 



where srzr(^) is the RZR contribution averaged over 
Rabi frequencies available in a MOT, as discussed in 



Sec. IIV Al The results of theoretical simulations accord- 
ing to formula ifTTHl are depicted in Fig. (fit 1°) and 
compared with experimental data. While the agreement 
of the modelled and experimental spectra improved con- 
siderably after considering both contributions to s{6), 
the theoretical spectrum still exhibits resonances that are 
too narrow. SimpHstic attempts to reduce resolution of 
the theoretical spectra by extending the fl inhomogene- 
ity range in the Srzr(<5) contribution does not improve 
the fit. Whereas the net spectrum s{6) becomes broader, 
its individual peaks and dips are shifted out of coinci- 
dence with the experimental ones. This indicates that 
some other kind of inhomogeneous broadening mecha- 
nism must affect the spectra. This mechanism is due to 
the quadrupole MOT magnetic field B and finite size of 
the atomic cloud. As the spatial distribution of atomic 
density in the trap is Gaussian, modelling of the indi- 
vidual RZR resonances, such as in Eq. Q, should be 
performed with the Gaussian profiles G{S, 6kj,B), rather 
than with the Lorentz functions iQ. We thus take 



2al 



with 



c^B = — ; ttCTz > lk,j- 

n oz 



(11) 



(12) 



In the above formulas, C is a constant, gp is the Lande 
factor of the ^Si/2(^' = 3) hyperfine state, hb is the Bohr 
magneton and ctz is the Gaussian radius of the atomic 
cloud along z axis. The results of the theoretical sim- 
ulation of the absorption spectra with inclusion of the 
MOT magnetic field are presented in Fig. as fit 2° and 
are compared with the experimental data. In contrast to 
averaging exclusively over |E| (fits 0° and 1°), the agree- 
ment of the theoretical and experimental absorption sig- 
nals is now very good. This indicates the importance 
of the inhomogeneous broadening of the MOT spectra. 
So far, the RZRs in a MOT have been considered to be 
free of any inhomogeneous broadening. Our experiment 
shows that this is not the case. While Doppler broaden- 
ing is negligible, the magnetic field inhomogeneities cause 
substantial broadening of the pump-probe Raman spec- 
tra and corresponding change of their Hneshapes from 
Lorentzian to Gaussian. 



D. Four- wave mixing signals 

In our experiment, the four- wave mixing spectra are 
recorded simultaneously with the probe beam absorption, 
which assures that they are acquired in exactly the same 
experimental conditions. The four- wave mixing signals 
recorded in our experiment (Figs. [JJd, lib) clearly exhibit 
two contributions: the broad pedestal and pronounced, 
ultra-narrow features for (5 w 0. The former can be at- 
tributed to the superposition of the RZRs and RIRs and 
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Figure 7: Experimental and theoretical absorption spectra for 
cr"*" and cr" -polarized probe beam. The experimental condi- 
tions: /ma>c = 10.2 mW/cm^ A = -SF, dBjdz = 13 G/cm. 
Fit 0° includes only inhomogeneity of the E-field, fit 1° - as 
fit 0° + RIR, fit 2° - as fit 1° + inhomogeneity of the B- 
field. Parameters of theoretical simulations: Rabi frequency 
averaging range € [SF; 7.5r]; Gaussian width of a single, 
inhomogeneously broadened resonance (Jb = O.OSSF, which 
corresponds to the actual trap size — 0.9 mm; the tem- 
perature of atomic cloud T = 155 /iK was fitted to achieve 
the best agreement of the simulation with experimental data 
and is consistent with the results of Ref. 23||. For better com- 
parison, the experimental spectra are added in grey to each 
fit. 



the latter is due to the transitions between kinetic mo- 
mentum states, RIR and possibly to the Rayleigh scat- 
tering Q. Ultra-narrow resonances are here far better 
resolved than in the case of corresponding absorption sig- 
nals. Hence, the four-wave mixing spectroscopy appears 
to be more sensitive to the recoil effects. However, precise 
verification of this conjecture is hampered by the level of 
complexity of the theory of four-wave mixing for recoil 
induced resonances (see, e.g. jl^l). Also, the nature of 
the four-wave mixing process differs substantially from 
the probe absorption: the four wave mixing signal is, in 
general, calculated as |pe3(p,p')P; where peg{p,p') is an 
off-diagonal element of a momentum-dependent density 
matrix. Since peg{p,p') = Peg{p,p')RZR + Peg{p,p')Rm, 
the interference terms Pegip,p')RZR x Pge{p',p)RiR play 
an important role in the four-wave mixing signal. These 
terms are likely to be responsible for significant differ- 
ences between four-wave mixing signals for the two cir- 
cular pro be polarizations. The theoretical work by Guo 
et al. |2Q| presents a derivation of the four-wave mixing 
resonances due to atomic recoil in a two-level system. 
However, for this theory to be appHcable to our case, it 
should be complemented with the RZR contribution for 
systems with nonzero angular momenta. 
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E. Optical anisotropy of cold atoms 

In Sec. IIV Al we discussed the RZR contribution to 
the MOT absorption spectra. It does not depend on 
the probe beam polarization, due to the fact that in the 
considered case of 7r-polarized pump field, neither popu- 
lations of individual Zeeman sublevels, nor the Raman- 
resonance frequencies, depend on the sign of mg. On 
the other hand, in the case of RIR, the polarization of 
the probe beam plays a crucial role. As was mentioned 
in the conclusions of Sec. IIV Bl recoil processes for two 
opposite circular polarizations of the probe result in dra- 
matically different momentum and energy transfers. The 
momentum transfer from trapping beams perpendicular 
to the probe is the same for both probe beam polariza- 
tions and, according to Eq. equals Ap — ±hk\/2/2. 
However, the probe beam polarization selects which of 
the two pump beams, the co-propagating or the counter- 
propagating with the probe, is involved in the RIR. The 
quantity of the momentum transfer is Ap « ±hk9 for the 
a~- and Ap « ±2hk for the cr """-polarized probe. Hence, 
the RIR contribution is the cause of the obvious difference 
between the absorption spectra for the two probe beam 
polarizations around (5 ~ 0. Fig. |Ht exemplifies such a 
difference, obtained by subtraction of two independently 
measured absorption signals, compared with the calcu- 
lated one (Fig. ISt). While wide structures of the two 
curves (between -0.1 F and +0.1 F) are well correlated, 
the agreement of their central parts is less satisfactory. 
The exact cause of this discrepancy is not yet well un- 
derstood. There are several possible effects that can in- 
fiuence that challenging recording of kHz-wide spectral 
features. On one hand the accuracy with which we de- 
termine the zero of the pump-probe detuning is limited, 
the relative phase of the two laser beams undergoes some 
residual fiuctuations and the MOT parameters may vary 
between successive acquisitions of the a~ and ct"'" spectra. 
On the other hand, additionally to the RZR and RIR con- 
tributions that we consider, the central structure might 
be systematically affected by additional mechanism such 
as, e.g. the Rayleigh scattering Q or Faraday rotation 
due to imperfect compensation of a magnetic field and / or 
imperfect balancing of the MOT beams. 

The asymmetry between absorption spectra, i.e. be- 
tween optical properties of the sample as seen by probe 
beams of different polarizations, is known as dichroism, 
or generally as optical anisotropy. In our case, the atomic 
sample is placed in a center of a MOT, where B = 0, 
and the optical pumping by pump beams creates atomic 
alignment, rather than polarization. Thus the optical 
anisotropy of our cold sample can be associated exclu- 
sively with the recoil processes. Consequently, when two 
absorption spectra for the cr~ and probe beams are 
subtracted, the RZR contribution is completely elimi- 
nated and one is left only with the contribution from 
RIRs for the pump beams co- and counter-propagating 
with the probe. 

The optical anisotropy associated with RIR and shown 
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Figure 8: The difference a~ —a for the experimental absorp- 
tion spectra presented in Fig.EI(a) compared with theoretical 
prediction (b). 

above (Fig. |SJ constitutes experimental evidence of the 
effect recently predicted by Dubetsky and Herman [^^l- 
The authors of this paper concentrated on the some- 
what different experimental situation in which the opti- 
cal anisotropy is detected via rotation of the polarization 
plane of a weak, linearly polarized probe beam propa- 
gating perpendicularly to the pump beam [^Jl- Still, the 
anisotropy described here has the same physical origin, 
the RIR. 



V. CONCLUSIONS 

We have successfully explained absorption spectra 
recorded in a working MOT. We have shown that full 
agreement between the theory and experimental data is 
attainable only when realistic MOT conditions, the light 
modulation in a trap and the quadrupole MOT magnetic 
field are taken into account. Our theoretical analysis em- 
phasizes the importance of the atomic recoil processes 
in a MOT. Recoil-induced resonances have already been 
used in atomic velocimetry [13J,[14J. In order to apply 
this method to atoms in a working MOT, one has to 
eliminate other processes infiuencing the spectra, namely 
the Raman transitions between light-shifted Zeeman sub- 
levels of the ground atomic state, RZR. For the MOT ve- 
locimetry based on RIRs, the RZR contribution appears 
as an undesired background which hinders precise deter- 
mination of the RIR width. Unlike the RIR contribution, 
that due to the RZR is subject to significant inhomoge- 
neous broadening by trap magnetic and light fields, so it 
is worthwhile to eliminate this background, e.g. such as 
suggested above by recording the optical anisotropy or 
the differential la~ absorption. 

Let us note that in contrast to the experiment by 
Schadwinkel et al. 0| with carefully phase-stabilized 
trap beams, we see no contribution from the Raman tran- 
sitions between vibrational levels associated with a pos- 
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sible optical lattice, whose frequencies could be in the 
similar range of the order of 100 kHz. We have found 
that under regular conditions in a standard MOT with 
nonstabilized trap beams, the short-term, random fluctu- 
ations of the trapping beam phases caused by mechanical 
instabilities of the setup, wash out the interference pat- 
tern, thus preventing stationary modulation of optical 
potential and atomic localization. 

Another aspect of our results is the four-wave mixing 
spectroscopy. Since the four- wave mixing signals are gen- 
erated in a more complex process, especially when the 
multilevel structure and recoil effects have to be taken 
into account, their theoretical analysis is complicated and 
requires further investigation. At the same time they ap- 
pear to be far more sensitive to recoil processes than ab- 
sorption, as the ultra-narrow structures are here highly 
pronounced. This interesting feature motivates one to 



work out the method of MOT diagnostics based on four- 
wave mixing. 
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